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CALCULATING- THRUST DISTRIBUTION AND EFFICIENCY 

OF AIR PROPELLERS . * 
By Theodor Bienen. 

Supplementary to various articles published by Professor 
Von Kan-nan and myself on the propeller theory,** I am now pro- 
posing a method for the preliminary approximate calculation of 
the thrust distribution and efficiency of air propellers under 
any operating conditions. 

The calculation is based on the following manner of rea- 
soning. The propeller is imagined to be replaced by a suitable 
mechanism, with whose aid a certain mass of air is driven oppo- 
site to the direction of flight (and in the peripheral direc- 
tion). Such a mechanism would, for example, be represented by 
a series cf blade elements which are warped with reference to 
one another in the radial direction from element to element, 
corresponding to the pitch and which, in the peripheral direc- 
tion, fill the propeller circle with infinite closeness. Every 

blade element generates a certain thrust, the magnitude and di- 

*"Die rechnerische Ermittlung der S chub vert eilung und des Wirk- 
ungsgrades fur ausgefuhrte Luf tschrauben bei beliebigen Betriebs- 
zustanden. " From "Zeitschrift fur Flugtechnik und Motorluft- 
schiffahrt, " Nov. 27, 1926, pp. 485-487.' 

**Bienen and Kan-nan, "Zur Theorie der Luf tschrauben, " Z.d.V.D-It 
("Zeitschrift des Vereines deutscher Ingenieure" ) , Nov. 29, and 
Dec. 20, 1924. Th. Bienen, "Die gunstigste Schubvert eilung fur 
die Luftschraube usw., 11 Z.F.M., May 28, and June 13, 1925. 
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rection of which are known (provided the performance coeffici- 
ents of the elements are known) and whose axial component rep- 
resents the contribution of that element to the propeller thrust. 

Hence our task becomes the following, namely, to determine 
the speed relations and the forces developed on a section with a 
given direction and Velocity of the air current, a problem which 
can be approximately solved with the aid of the momentum theory. 

There still remains to be considered the effect of the fi- 
nite number of blades. This is determined in the following man- 
ner, as suggested by Professor Prandtl. 

Effect of the Number of Blades 

In accordance with our assumption, the thrust at any defi- 
nite distance from the center of the propeller circle is uni- 
formly distributed throughout the circumference, while the thrust 
distribution along the radius can be assumed to be distributed 
as raa,j be desired. 

Figure 1 shows the speed relations of the propeller blade 
at the distance r from the propeller axis. 

w = flight speed, 

oo = angular velocity of propeller, 

= arc tan the effective angle of pitch, 

w +"!' ' 

(3 r ' = arc tan £ — tc, the induced angle of pitch. 

rf a) - Jf\ 

c' = the added (or induced) velocity. Its components in 
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the axial and peripheral direction are respectively 
w' and rw'. (The induced velocity first reaches 
half its value at the position of the propeller 
"blade itself . ) 

With a finite number of blades, there is a decrease in the 
thrust at the blade tips, which is a natural result of the re- 
lease of an especially large number of vortices at those points. 
The velocity distribution in the periphery is no longer uniform. 
If the periphery is developed and the induced velocity is plotted 
over it, it is then found that it attains its maximum value at 
the propeller blade, but falls to a certain value before and be- 
hind the blade. 

The mean induced velocity was approximately calculated by 
Prandtl* as w mean = k w 1 , in which w 1 denotes the maximum 
effective velocity at the propeller blade and k equals 

2 arc cos e _ R^jr ; 

in which z denotes the number of blades and (3' the induced 
angle of attack at the blade tip. The propeller with a finite 
number of blades can now be replaced by a dynamically equivalent 
propeller with an infinite number of blades. Prandtl computes 
the reduced radius of this propeller by the formula 

H red . = R (1 - 1^^*} 

*See A. Betz, "Schraubenpropeller mit geringstem Energieverlust" 
with Appendix by L. Prandtl, "Nachr. d. Ges. d. Wissensch. zu 
Gottingen," "Math.-physik. Kl." 1919, p. 193 ff. The quantity k 
is given in the form of a curve in the previously mentioned paper 
by the writer. 
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A knowledge of the maximum velocity at the propeller 
blade, as also of the mean velocity, is very important for de- 
termining the shape of the propeller. Prandtl' s correction 
factors seem to approximate the actual relations closely enough. 
At least the agreement was excellent between them and the ac- 
tual velocity and. thrust distribution obtained in the wind- 
tunnel experiments at the Aachen Aerodynamic Institute. 

Calculation of the True Angle of Attack 

As shown in Fig. 1, the task is only to determine the dif- 
ference between the geometrical and the induced angle of attack. 
The induced velocity c' , with its axial and peripheral compo- 
nents, can be calculated directly from its geometrical rela- 
tions. With the induced velocity we then know the thrust dis- 
tribution and the efficiency. 

If we apply the momentum theory to an annular area, of 
width dr, in the propeller circle at the distance r from 
the propeller axis, the thrust developed by the propeller ele- 
ments, perpendicular to the direction of attack, can be written 
as follows: 

mo' * = A 

1 f w + |') 2 r tt dr c» k = | ^ c a z t dr (l) 

in which c'/2 is the maximum value of the induced velocity; 
k, the Prandtl correction; c a , the lift coefficient of the 
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"blade section; and zt, the sum of the blade chords of all the 
sections comprised in the radius under consideration. The re- 
maining quantities, in so far as they are not given in Figure 1, 
are designated in the usual way. 

On the basis of the geometrical relations in Figure 1, we 
obtain 

c« = 2 v sin 6 = § w n sin 6 

sin p 

v' = v cos 6 = * - cos 6 

^ (2) 

w« = c' cos 0» = 2 . W Q sin 6 cos (P + 6) 

sm p 

r 0)' = c' sin P' = ? w n sin 6 sin (P + 6) 

sm p 

If i denotes the angle between the line of zero thrust 
( ca = 0) of the section and the "effective" direction of attack 
(i.e., of the apparent direction given by the operating condi- 
tions), we can write, in the usual manner, 

c a = ca' (i - 6) (3) 

in which ca' is derived from Ca , according to the angle of 
attack. In what follows, we will treat c a ' as a constant. 
For ordinary blade sections and the usual angle-of-attack re- 
gion, this approximation answers very well. We can then put 
c a ' = 5.25 (the angle in circular measure). In general, 
c a ' has smaller values only at very great angles of attack. 
This should be remembered when the following conclusions are 
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applied to the operation of propellers under conditions in which 

the above assumption is not (or is only partially) fulfilled. 

test 

This case will occur especially with the propeller on the/ stand 
and there more often with propellers having a large V/NO ratio 
(i.e. , large ratio of flight speed to peripheral velocity, hence 
in low-speed propellers). Since, however, this disadvantage is 
chiefly confined to the elements near the hub, the resulting 
error is small, for these blade elements furnish only a small 
fraction of the total thrust. 

From equations (l), (2), and (3), we now obtain 

f -, , sin 6 , Q , „ v ^ A . z t cos 3 6c a '(i- 6) 

Since 6 is always very small, we expend sin 6 and cos 6 . 
By neglecting everything after the first term, we obtain the 
difference between the induced and effective angle of pitch 
(angle of slip) . 

i 



1 + K 8 t n J sin P, 



(4) 



c a ' zt 

Taking into account the quadratic .terms .of 6 we .obtain 



K Birr sinpj + ^ KStTr cos|3 r + ^ 



K 



8rrr c 



2L§c + i . ^sirr B mp r + ^ 



z t c a ' >J V z t c a ' / 



Even this formula, with a suitable tabular arrangement, 
can be evaluated much simpler and quicker than perhaps appears 
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possible at first glance. A n y one desiring greater accuracy- 
can make a second calculation with the more exact value of c a ' 
He will then obtain accurate values even for points at any de- 
sired distance from normal operating conditions. 

Calculation of Thrust Distribution and Efficiency 

The induced velocities in the axial and peripheral direc- 
tions can be calculated with the aid of equation (2). From the 
momentum and blade-area theories, in the usual manner for the 
annular element with the radius r and the width dr, we then 
obtain the contribution to the thrust 

<|S = 7 / w + w_V(l- £ tg p r i) 2 hit 
dr g V 2 / to x 

and the contribution to the torque (6) 

M = 1 { * + ™ l \ r = w « (i + € ctg P r O 2 r tt 
dr g v. 2 / or' 

in which € = ^3E, the D/L ratio Of the element (for an inf i- 
°a 

nite aspect ratio). Since the angle of attack is given by 6, 
both c a and c v/ are known for the sections tested as models. 
The D/L ratio can be approximately determined for elements 
for which no results of model tests are available. We will re- 
cur to this point later. 

The thrust and torque distribution along the radius is 
given by equation (6). The total thrust and torque are ob- 
tained by integration. 



N.A.G.A. Technical Memorandum No. 444 



8 



The efficiency, at the distance r from the propeller axisj 

is 

- = d S w (n\ 
We finally obtain for the total efficiency 

/ *r 

This integration can be quite easily made, either numerically 
or graphically. 

Remarks on the Method of Calculation 

The calculation is not difficult when test results of the 
individual propeller elements are available, which is rarely 
the case- For calculating 6 from equations (4) and (5), we 
must know the c a = 0 line (or the angle i) and the c a and 
c w values in terms of the angle of attack (i - 6). For the 
approximate determination of the c a = 0 line of a known ele- 
ment, the following method is more suitable.* 

Connect the trailing edge of the elements (Fig. 2) with the 
leading edge (contact points of a circumscribed circle) by the 
straight line a. Next draw the line b, so as to bisect the 
angle formed by the tangents ''o the sides of the trailing 

edge.. The straight line bisecting the angle between a and b 

*See Karman and Trefftz, "Potentialstromung um gegebene Trag- 
f lachenquerschnitte, " Z.F.M. 1918, p. Ill ff. 
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is the c a = 0 line (or line of zero thrust). This is exactly 
true only for elements having. a trailing-edge angle of 0, 
i.e*, for so-called "Joukowsky profiles." For finite trailing- 
edge angles, it would be necessary to add a correcting term de- 
pendent on the magnitude of this angle. This term is, however, 
so small as to "be practically negligible. 

If the same basic profile is employed for the different 
cross sections by geometrically distorting it, then a single 
determination of the c a = 0 line is sufficient. For any other 
thickness ratio d/t, the angle between the c a = 0 line and 
the chord is 

a - 4 

in which a© and f^'^ are the corresponding values, either 
measured or otherwise determined. For any desired angle of at- 
tack c a is then found according to equation (3). 

In determining the c w values, care must be exercised that 
they represent only the coefficient of the pure profile drag 
(i.e., without induced drag). The profile drag may, however, 
with sufficient accuracy, be regarded as constant in the custom- 
ary angle- of-attack region. It can be estimated by comparison 
with known profiles. . 
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Summary 

A method is presented whereby the thrust and torque dis- 
tribution, the efficiency, and therewith the requisite engiri'e 
power, can be determined for any propeller under any working 
conditions. With the aid of the momentum and blade-area theo- 
ries., on the one hand, and the results of the elementary air- 
foil theory, on the other hand, the only unknown quantity, 
namely, the actual aerodynamic angle of attack of each propel- 
ler-blade element is determined. Herewith all other quantities 
can be calculated in a simple manner. 

Translation by Dwight !!• Miner, 
National Advisory Conmitt 
for A eronautics. 
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Figs. 1,2. 




Line of c a = 0 



(zero thrust) 



a 



f — r: ^'^^P^^ 



Angle between tangents 
at trailing edge. 



Fig. 2 .Graphic determination of line of zero thrust, 
(c a =0). . 



